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We study macroscopic observables defined as the total value of a physical quantity over a collection 
of quantum systems. We show that previous results obtained for infinite ensemble of identically 
prepared systems lead to incorrect conclusions for finite ensembles. In particular, exact measurement 
of a macroscopic observable significantly disturbs the state of any finite ensemble. However, we show 
how this disturbance can be made arbitrarily small when the measurement are of finite accuracy. 
We demonstrate a tradeoff between state disturbance and measurement coarseness as a function of 
the size of the ensemble. Using this tradeoff, we show that the histories generated by any sequence 
of finite accuracy macroscopic measurements always generate a consistent family in the absence of 
large scale entanglement, for sufficiently large ensembles. Hence, macroscopic observables behave 
"classically" provided that their accuracy is coarser than the quantum correlation length-scale of the 
system. The role of these observable is also discussed in the context of NMR quantum information 
processing and bulk ensemble quantum state tomography. 

PACS numbers: 03.65.Ta, 03.67.-a 
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I. OVERVIEW 



Macroscopic observables correspond to physical quan- 
tities which are accessible to our senses. Since the phys- 
ical scale of individual quanta is generally tiny, macro- 
scopic observables arise when a collection of quantum 
systems are measured jointly. Formally, they can be 
described by type projectors, which reveal information 
about the average population of single-particle states. 
For example, the total magnetization of an ensemble of 
spin-i particles provides some information about the rel- 
ative occupation number of the up and down states. We 
will derive several general properties of these measure- 
ment and discuss how they lead to the emergence of a 
quasiclassical domain in the absence of large scale entan- 
glement. 

The effect of macroscopic observations on infinite en- 
semble of identically prepared quantum systems has been 
studied in various contexts 0, 0, El El • The main conclu- 
sion of these studies is that the state \vj)® N describing 
such an ensemble is an eigenstate of type projectors when 
N = oo. However, for finite ensembles, things change 
dramatically. The measurement of a macroscopic observ- 
able induces a disturbance which increases as the size 
of the ensemble grows, in apparent contradiction with 
the infinite-copy result. This discrepancy follows from 
the ambiguous extension of finite-copy considerations to 
the nonseparable Hilbert space of an infinite-copy ensem- 
ble In this article, we show how the essence of the 
infinite-copy result can be recovered by "smoothing" the 
type projectors into coarse-grained positive operator val- 
ued measurement (POVM) (essentially going from the 
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strong to the weak law of large numbers) . 

The paper is organized as follows. The central math- 
ematical objects of the present study are defined in Sec- 
tion^ We first summarize the method of type and define 
type projectors. These are projectors on the degenerated 
eigensubspaces of macroscopic observables of the form 
An = X^fcLi a {k) i where ar^ acts on the kth system of the 
ensemble. Using the theory of generalized measurements, 
we also define coarse-grained POVMs corresponding to 
finite accuracy estimation of a macroscopic observable. 

Section IfTTI contains the core mathematical analysis of 
our study. We first recapitulate the well known facts 
about type projectors acting on infinite ensemble and 
show how they dramatically break down for finite en- 
semble. Then, we show how the result is approximately 
recovered when the measurements are of finite accuracy, 
and study the general tradeoff between measurement 
coarseness and state disturbance — measured in terms 
of fidelity — as a function of the size of the ensemble. In 
short, we demonstrate that a measurement of coarseness 
(7 l/v^/V leaves the systems essentially unchanged, i.e. 
the fidelity F between the pre- and post-measurement 

state of the ensemble satisfies 1 — F oc ln ^.2 - ■ 

Section El is a discussion of the de Finetti representa- 
tion theorem which provides a wide class of state — ex- 
changeable states — for the study of macroscopic observ- 
ables. Exchangeable states have recently been employed 
for the discussion of quantum state tomography based 
on single-system measurements followed by Bayesian up- 
date [6(. We will show how macroscopic observables of- 
fer an alternative perspective on quantum tomography. 
Moreover, this approach offers interesting applications 
for quantum information theory 0, an( I is a more accu- 
rate description of experimental spectroscopy-based im- 
plementations of tomography, e.g. as achieved in Ref. 
Macroscopic observables also provide an explanation 



2 



for the emergence of the classical world we perceive from 
the underlying quantum theory. Indeed, we demonstrate 
in section E] that in the absence of large-scale entan- 
glement, one of the main characteristic of the classical 
domain — namely the commutativity of classical ob- 
servables — follows naturally from the general proper- 
ties of coarse grained type POVMs. This is done us- 
ing the consistent histories formalism which we briefly 
summarize. We demonstrate that the histories generated 
by any sequence of macroscopic observables of accuracy 
u y£jN are consistent, where £ is the quantum corre- 
lation length-scale of the system. This generalizes some 
ideas introduced by Halliwell 10] on how to achieve clas- 
sicality in closed quantum systems. 

Finally, Section IVII discusses the role of macroscopic 
observables in NMR quantum information processing. In 
this context, macroscopic observables are used to extract 
the output of the computation, but also, since the mea- 
surement device can not be "turned off" — i.e. the state 
of the processor can always be read-off from the spec- 
trometer — they constantly perturb the computation. 
Following the results of Section IIIII and a measurement 
model introduced in Ref. , we show that the measure- 
ments used in NMR can in principle be sufficiently pre- 
cise to extract useful information about the computation 
but yet so coarse grained that they induce a negligible 
perturbation. However, as we will demonstrate, NMR 
measurements may not follow our optimal measurement 
coarseness-state disturbance tradeoff when performed at 
room temperature; caution is advised when applying our 
conclusions. 



II. DEFINITIONS 

This section contains all the mathematical definitions 
required for our study. For sake of clarity, we adopt 
the vocabulary of NMR. Our general setting consists of 
an ensemble of N quantum systems of the same nature. 
Therefore, we shall refer to individual systems of an en- 
semble as molecules and to the ensemble of N molecules 
itself as the sample. Thus, the word "molecule" should 
not be taken literally in what follows; it could be any 
elementary constituent of a larger system. 



A. Method of types 

The method of type is a very powerful statistical tool 
with applications ranging from large deviation theory, 
universal coding, and hypothesis testing. We will only 
scratch the surface of this theory here, more details and 
applications can be found in [l2j for instance. 

Let X = Xj ± Xj 2 . . . Xj N £ X N be a string of N let- 
ters drawn from a d-letter alphabet X — {xi, x 2 , ■ ■ ■ Xd}- 
The type (or empirical probability distribution) of X is a 



vector of positive numbers summing to one defined by 

L{X)=(L 1 (X),L i {X),...,L d {xj), 

where Lj(X) is the relative frequency of the letter Xj in 
the string X 

1 N 

fe=l 

it is simply the number of occurrences of the letter Xj 
in X, divided by the length of X. For example, if X — 
{a,b,c} and N = 4, then L(cbaa) = (^, \, \). We also 
define a type class T to be the set of strings of a given 
type: 

T[L] = {X eX N : L(X) = L}. 

For example, using the same alphabet as above, we have 
T[(i 0, 1)] = {accc, cacc, ccac, ccca} . The class T[L] can 
be generated by applying all permutations to any single 
string of type L. Hence, the number of elements in T[L] 
is given by the multinomial coefficient 

|T[L]| =( N ) 

1 L Jl \NL 1 ,NL 2 ,...,NL d J 

Nl 

~~ {NL 1 )\,(NL 2 y.,...,(NL d )\ 

Let R = (Ri, i?2, ■ ■ • , Rd) be a probability distribution 
over X. The probability of the string of outputs X = 
Xj 1 . . . xj N of N letters, each drawn independently ac- 
cording to the distribution R, is P(X) = Rj 1 Rj 2 . . . Rj N . 
This can also be written as 

F(I),C m < i2m ...C (I », 

so given a fixed distribution R, the probability of a string 
X 6 X depends only on its type. Intuitively, the type 
of the observed outcome X is very likely to be close to 
the probability distribution of the random variable, i.e. 
Lj(X) ?v Rj : as N increases. This is the substance of the 
typical sequence theorem 

P(||L - R||? > e) < ~ e -Ne/2 (1) 

where the "difference" between the type L and the prob- 
ability distribution R is quantified by the variational dis- 
tance (Li-norm) 

\\L-B.\\ 1 = J2\L j (X)-R j \. 

3 

The typical sequence theorem takes on various forms. It 
can be formulated in a stronger version using the relative 
entropy, which is an upper bound to the variational dis- 
tance. Nevertheless, for our considerations, this simple 
version will be sufficient. 
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B. Macroscopic observable 

Using this notation, we now formally define macro- 
scopic observables. Consider a Hermitian operator (i.e. 
observable) a acting on the d-dimensional Hilbert space 
of a single molecule TL m - Let {|xi), \x2), ■ ■ ■ \xd)} and 
«2) ■ ■ ■ j a d} denote its eigenvectors and eigenvalues: 
a\xj) — a.j\xj). We will assume that a is non- degenerate, 
generalization is straightforward. The macroscopic ob- 
servable An corresponds to the sum of observable a over 
all the N molecules of the sample, 



and Q^!'^ associated to two distinct macroscopic observ- 
ables An = ^ k anA and Bn = J2k fyfc) n °t commute, 
unless the underlying single-molecule observables a and b 
happen to commute. To avoid cumbersome notation, we 
will only use an extra superscript when necessary (c.f. 
Section llV Afl . For the time being, we will consider a 
fixed arbitrary macroscopic observable An- In this case, 
it is straightforward to verify that the type projectors are 
mutually orthogonal and that they sum to the identity 



Qi N) Qi N) 



11. 



(4) 



N 



An — ^ ei(fc), 



(2) 



k=l 



where d(u\ is the operator a acting on the fcth molecule, 
O(fe) = H ® . . . ® 11® a ® H ® . . . (g) 11. 



fc-i 



N-k 



The operator An acts on the joint Hilbert space of the 
N molecules Ti s = TC^ N — the Hilbert space of the sam- 
ple — which has dimension d . We use the standard 
abbreviation \X) = {x^) <® \xj 2 ) ® • • • <8 \xj N ) for each 
string X € X . Clearly, the states \X) form an or- 
thonormal basis for Ti s . Moreover, they are eigenstates 
of the macroscopic observable An- 

A n \X) = A N \x h ) ® \x h ) ® ...0 \x jN ) 

]a (k) \x 3l }(g>\x j2 }® \x jN ) 



N 



fc=l 
N 

fe=l 



ajjxjj ® \x h ) ® . . . ® \x jN ) 



N 

£«*» 

k=l 



\X)= [^NLjiXfa ) \X) 



Thus, we see that the eigenvalue associated to a basis 
state \X) depends only on its type L(X). As a conse- 
quence, the degenerated eigensubspaces of An are those 
subspaces spanned by the vectors \X) belonging to the 
same type class. 

This brings us to the definition of type measurements 
which are von Neumann measurement composed of the 
projection operators on the subspaces of a given type: 



QW= £ \X){X\. 



(3) 



XGTVL] 



Each of these type projectors is labeled by a vector of d 
positive numbers Lj which correspond to the type L(X) 
of the basis vectors \X) spanning the subspace. Obvi- 
ously, the projectors depend on the choice of basis 
\xj) over Ti m , i.e. on the eigenvectors of the observ- 
able a, so we could explicitly note Q^ ,a \ Moreover, we 
would like to stress that the spectral projectors Q^ ,a ^ 



In words, these projectors correspond to the exact mea- 
surement of the population of the levels \xj) over an en- 
semble of N molecules, without distinguishing between 
the molecules of the sample. The type projectors Qjf^ 
allows us to express the operator An in a simple form: 



A 



N 



(5) 



where we have defined A-^, — 53j=i NLjOtj. Similarly, 
any macroscopic observable of the form Eq. @ has a spec- 
tral decomposition involving only type projectors, as in 
Eq.©. Hence, following textbook quantum mechanics, 
when measuring a macroscopic observable — or mea- 
suring the "expectation value" of a physical observable 
over a macroscopic sample — , one is really performing a 
projective von Neumann measurement composed of type 
projectors. 

These type projectors have been studied under many 
different forms and take on many different names. 

Among other formulation are the frequency operators. 
Recall that Lj(X) is the relative frequency of the symbol 
Xj in the string X . We can define a frequency operator 

FjW = <£ L .( X )\X){X\. 
x 

This operator is a macroscopic physical observable whose 
eigenvalues are fj = 0,-^,-^,...,l. Indeed, iv takes 
on the form of Eq. (|2J) by setting the single-molecule ob- 
servable a to jj\xj)(xj\. Following textbook quantum 

mechanics, when the measurement associated to Fj N ^ is 
performed and eigenvalue fj is observed, the state of the 
system gets collapsed to the subspace spanned by the 
states \X) for which Lj(X) = fj. Hence, the eigenvalue 
fj indicates the relative population of the single-molecule 
state \xj) in the sample of A" molecules. 

The above construction yields d commuting physical 
observable {Fj }j=i,...,d, one for each single-molecule 
state {\xj)}j=i t ,.. t d- Regrouping these observable into a 
d-component observable yields 



P(N) 



(AD p (JV) 



, F, 



(6) 



(AT) 

L ) 
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which takes on the form of Eq.JSJ, with a d-component 
eigenvalue = L. The value of any macroscopic ob- 
servable of the form Eq.@ can be deduced straightfor- 
wardly from the value of F^. Hence, a great deal of 
attention has been focused on the macroscopic observ- 
able without loss of generality. 

We illustrate macroscopic observables for a sample of 
N spin-i particle. We choose the basis \x%) = |"f) and 
I £2) = II) corresponding, respectively, to +^ and — i 
units of magnetization in the z direction: 



where the Kraus operators Aji can be any set of oper- 
ators satisfying Yli^-ji-^ji = Ej- Here, we will often 
consider ideal quantum measurements where the distur- 
bance inflicted to the system is minimal. This restriction 
is necessary if we want to study the optimal tradeoff be- 
tween information gathering and state disturbance. To 
each measurement outcome E, of an ideal measurement 

Ej . In this 



is associated a single Kraus operator Ajo = 
case, the state update rule Eq.© simplifies to 



<j z \xi) = and a z \x 2 ) = -^\x 2 ). 

We can use a single positive number L € {0, ■h, j|, . . . , 1} 
to label a type of a binary string X, which corresponds 
to the fraction of Xi's (or spin up's) in X. Hence, a type 
L is a shorthand for L = {L, 1 — L). The bulk (or to- 
tal) magnetization of the sample is equal to the sum of 
the magnetization of each molecules: the corresponding 
operator is therefore M z = J2k=i a \ky wnere a? k -. is the 
"z" Pauli operator acting on the fcth molecule. When the 
sample is in a state of a definitive type, its bulk magne- 
tization is equal to \X(L\ — L 2 ) = \N(2L — 1), which is 
simply the number of spins pointing up minus the number 
of spins pointing down, times ^- Hence, the observable 
corresponding to the bulk magnetization can be written 
as 



L 



where the sum is over all types. The type projectors 
Q^iP are projectors on the degenerated eigensubspaces 
of the bulk magnetization operator. Clearly, an exact 
measurement of the magnetization M z would reveal the 
type of the state of the sample, i.e. the relative frequency 
of up and down spins. 



C. Coarse grained macroscopic POVMs 

We will now present how finite accuracy macroscopic 
observables can be expressed in terms of type projectors. 
Before we do so, we briefly recall some basic concepts 
of the theory of generalized measurement. Generalized 
measurements (POVMs) are described by a set positive 
operators Ej summing to identity. The generalized Born 
rule for the probability of getting outcome Ej given initial 
state p is the same as for von Neumann measurements 



P(E J \p)=Tr{Ejp\. 



(J) 



After the measurement outcome Ej is observed, the state 
of the system gets updated to 



P(Ej\ P ) 



(9) 



which reduced to the regular state update rule when Ej 
are projection operators. Hence, von Neumann mea- 
surements are POVMs with an extra orthogonality con- 
straint. Any generalized measurement can be realized 
physically by coupling the system of interest to a larger 
system and performing a von Neumann measurement on 
the larger system; an example of such a physical con- 
struction will be presented in Section IVII Similarly, any 
such "indirect" measurement corresponds to a POVM. 
Hence, POVMs do not add anything extra to plain text- 
book quantum mechanics, beside conciseness. 

Continuing with our example, finite accuracy measure- 
ment of the bulk magnetization of a sample of N spin- 
i molecules can be described in terms of coarse-grained 

type operators ■ When the state of the sample is of a 
definite type L, the observed value of the bulk magnetiza- 
tion will not necessarily be equal to ^N(2L — 1) but, due 
to the uncertainty of the measurement apparatus, may 
take different values ^N(2£ — 1), with respective prob- 
abilities The function qi,{l) should be centered 
around L and have a certain width a corresponding to 
the coarseness of the measurement. 

Hence, the coarse-grained type measurements can be 
defined by "smoothing" the exact type projectors: 



(10) 



where qi_,(£) is some probability distribution over t cen- 
tered roughly at L and has the interpretation given 
above. In principle, £ could be any real d-dimensional 
vector, as it contains statistical fluctuations. For exam- 
ple, <tl(£) could be a (i-dimensional Gaussian 



1 



2n<T 2 



exp 



L 



2a 2 



(11) 



which is properly normalized J qj_,(£)d£ — 1 and where 
the L2-norm is \\£ — L||| = J2j(^j ~ Lj) 2 . The opera- 
tors Ei — Q^iQi^ form a POVM (with a continuous 
outcome) since they are all positive operators and satisfy 



Y.;^,P-\n 
P(Ej\p) 



(8) 



Q^Q^d£ = K. 



(12) 
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These coarse-grained type operators describe a situation 
where our measurement apparatus is not sufficiently pre- 
cise to measure the exact population of each level, but 
rather provides an estimation of it within a finite accu- 
racy a. 

We have assumed that the measurement outcome I 
takes on a continuous spectrum. However, several mea- 
surement apparatus, like those equipped with a numeri- 
cal output display, have a discrete spectrum of outcomes. 
This can be taken into account by choosing a smoothing 
function 



where {(.j} is the set of possible outcomes. Thus, we will 
henceforth consider the more general continuous case, 
but all our analysis carries through for discrete measure- 
ment outcomes by performing the above substitution. 



III. TYPE MEASUREMENT ON IDENTICALLY 
PREPARED SYSTEMS 

Type projectors were first studied by Finkelstein 0, 
Hartle Q, and Graham Q as part of discussions on 
the interpretation of probabilities in quantum theory. 
The main characteristic of type projectors identified by 
these authors can be summarized as follows. Let \tp) = 
S j Pj\ x j) be an arbitrary pure state of a ci-level molecule, 
with associated density matrix v = \ ?p |. Consider a 
sample of N identically prepared molecules, such that the 
state of the sample is I'fjv) = \ip)® N . Upon measurement 
of the type of the sample, we expect a result close to the 
probability distribution R = ((aJil^lxi), . . . , {xd\v\xd)) = 
(|/3i| 2 , . . . , |/3d| 2 )- Indeed, it follows from the strong law 
of large numbers that 



limit does not uniquely defined the operator on the non- 
separable Hilbert space Ti. m (g> TL m ® . . . of the infinite 
sample: specifying the action of F^°°^ on all states of the 
form \x%) <g) \x2)® ... is not enough to define it. This 
was realized in |j] where an alternative derivation of the 
probability rule was presented. Nevertheless, the pro- 
posed solution is still not satisfactory as it relies itself on 
probability theory. An up-to-date and rather critical dis- 
cussion of the status of the frequency operator and the 
related programs can be found in a recent paper of Caves 
and Schack [|| . We will not address these issues any fur- 
ther and do not claim to offer an alternative derivation 
of Born's rule. 

The second difficulty which is directly relevant to the 
present study concerns state disturbance. When a sys- 
tem is prepared in an eigenstate of a physical observ- 
able, the act of measurement does not disturb it. While 
Ea. (|f 3|) does not grant this for any finite N, one natu- 
rally expects that, as N grows, the disturbance caused by 
the measurement should decrease and eventually become 
negligible for all practical purpose. 

In what follows, we will show that the measurement 
of macroscopic observables induces an importance distur- 
bance to the state of the sample. In fact, this disturbance 
increases as the size N of the sample grows. This is in ap- 
parent contradiction with the conclusion that one might 
intuitively draw from Ea. l|13|) by extending it to finite N. 
However, we will show how the above conclusion can be 
recovered when the measurement of macroscopic observ- 
ables are of finite accuracy: sufficiently coarse grained 
type measurements induce a negligible disturbance to the 
state of the sample. We are interested in the tradeoff 
between measurement accuracy and state disturbance. 



A. State disturbance 



lim 

JV— >oo 



fW\* n ) 



\0i\ 2 \* N ) 



= 0, 



(13) 



where Fj is the jth component of the frequency opera- 
tor defined at Eq.®. In other word, F<W|#jv) =R|*at) 
with probability one in the limit of infinite N. This led 
Hartle to the conclusion that an infinite number of identi- 
cally prepared molecules are in an eigenstate \^ao) of the 
frequency operator F^°°' with eigenvalue R. Finkelstein, 
on the other hand, concluded from Eq. (| 1 31) that for finite 
N, \^n) is "close" to an eigenstate of F^' with eigen- 
value R. Thus, a measurement of the frequency operator 
reveals the probabilities Rj = (xj\v\xj), in the standard 
Copenhagen sense, of observing a single molecule of the 
sample in the state \xj). 

However, the conclusions reported above can be quite 
misleading. There are really two distinct issues here. 
The first one concerns the validity of the argument as a 
derivation of Born's rule to assign probabilities in quan- 
tum theory. The main complication comes from the def- 
inition of F(°°) as the limit of a finite operator. This 



The state of the sample can be rearranged as follows 



®JV 



I* 



5> 



e n* 

xex N \i=i 
d 



NLj(X) 



\X) 



E 



n/r' 



4=1 
d 



E m 

X£T[L] 



= E 11^4 I VW5l|L) (W) 
l \4=i 

where we have defined the normalized state 

il) = - 1 y \x) 
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and |T[L]| = ( A , L N NL ) denotes the cardinality of the 
type class T[L]. The density operator associated to this 
state will be denoted p N = |*jv)(*jv| = v® N . 

Upon measurement of the coarse-grained operators of 
Eq. (|10|l . the probability of observing an outcome within 
an infinitesimal volume range d£ of £ is P(Q[ N) \p N )d£ 
where (see Eq.Q) 

P(Q[ N) \ P n) = Tr{Q^p N Q { P) (15) 
= E V^M^W)Tr{Qi N) PN Qi N) } 

L,L' 



>N 



(16) 



and m(L,R) denotes the multinomial distribution 
(jVii N NL ) Ilj R^ Ll ■ Following Eq.©, the conditional 
post-measurement state of the ensemble given measure- 
ment outcome £ is 



PN\t 



eeii/t^ 



The post-measurement state is obtained by averaging the 
conditional post-measurement states over all measure- 
ments outcomes 



Pn = I P{Qf ] \pN)PN\idi 



Q { e N) \nm 
= EEIK L ^' 

L L' j,j> 

xG(L,L'W\T[L]\-\T\L']\ |L)(L'| (18) 



(A f ) 



where we have denned the decoh 



erence 



kernel 



G(L,L') = / y/ (£)di. 



(19) 



Notice that setting G(L,L') = 1 in Eq.fTHJ would yield 
a density matrix p' N identical to pn- Finally, the post- 
measurement state of a single molecule of the sample is 
obtained by taking a partial trace over N — 1 molecules 
p[ = TrN^i{p' N }, and similarly for the conditional post 
measurement state px\p = Trjv-i{/?jv|f}. 

The disturbance caused by the measurement is eval- 
uated with the fidelity between the pre- and post- 
measurement state. A fidelity of 1 indicates that the 
two states are identical — i.e. the measurement did not 



cause disturbance — while a fidelity indicates maximal 
disturbance. The fidelity between two states p and v is 



F{ P ,v) 



Tr 



p2 l/pl 



(20) 



If one of the state is pure, say v — \(j)){4>\, this reduces to 
the familiar "overlap" F{p,\<j)){4>\) = ((j>\p\<j>). 

It is instructive to first consider the case where the 
measurement are perfectly accurate, a = in Eg. (1111) . 
which implies qj_,(£) = S(£ — L) and G(L, L') = <5ll'- In 
this case, the post measurement density matrix is 



Pn 



£>(L,R) |L)(L| 



(21) 



it has completely decohered in the type basis |L), i.e. 
there are no off-diagonal terms of the form |L)(L'| like 
in Eq. I|18|) . The fidelity between the pre- and post- 
measurement state is then 

F„= (pn,p' n ) = ^[to(L,R)] 2 

L 

— m(L, R) x ^maxm(L,R)j , 

(22) 



(27riv)^rLi& 



where the subscript a = indicates that the measure- 
ment are perfectly accurate, and we have used Stir- 
ling's approximation in the last line. Clearly, exact type 
measurements greatly disturb the system, since fidelity 
goes to zero as the size of the sample increases, ex- 



cept in the case where /3 7 



J 3o ■ 



A similar conclu- 



sion based on different considerations was reached by 
Squires ^ follows from the concavity of fidelity 

F(P> J2j Vj v j) > Ylj PjF{p, Vj) that the conditional post- 
measurement state Pn\i has, with high probability, a van- 
ishing fidelity with the original state pn- 

The disturbance caused by an exact type measurement 
is most obvious when considering the conditional post- 
measurement state of a single molecule from the sample. 
As shown in the Appendix, 



Pi\i 



3 = 1 



(23) 



the conditional state of a single molecule is diagonal in 
the \xj) basis with eigenvalues given by the observed 
type of the sample £, independently of its state v prior 
to the measurement. However, following the typical 
sequence theorem Eq. (QJ, the observed coefficients tj 
are very likely to be close to (xj\v\xj). When averag- 
ing over measurement outcomes, we recover the state 
v — > P x\i = J2j Rj\ x j)( x j\ which has no off-diagonal 
terms, i.e. |afjX#,j|- Thus, the exact measurement of a 
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macroscopic observable completely decoheres individual 
molecules of the sample; it leaves the diagonal elements 
of v unchanged while suppressing all off-diagonal terms. 
(This situation might appear worrisome for bulk-state 
quantum computing; we will return to this in Section lVlJl . 
Moreover, the measurement creates correlation between 
the molecules, so p' N ^ {p'x)® N and p N \ t ^ (pi\e.)® N in 
general. 



B. Gaussian smoothing 

We now turn our attention to the case where the 
smoothing function qi,(£) has a finite width a. In the 
case of interest, the initial state of the sample \^>n) is 
pure, so combining Eqs. i|14|) and (|18|l we get 



F{pn,Pn) = {^n\p' n \^n) 



(24) 



L,L' 



m(L,|/J 3 | 2 )m(L',|ft| 2 )G(L,L') 



For sake of definiteness, we will consider the Gaussian 
distribution qj_,(£) defined at Ea. (|ll|l . The decoherence 
kernel defined at Eq. H19(l is then given by 



G(L,L') 



1 



27TCT 2 



-di 



exp 



2(2af 



This is not surprising as the decoherence kernel is the 
convolution of the smoothing function with itself. The 
convolution of two distribution of width o\ and 02 gives 
a distribution of width a = g\ + 02, so G(L, L') is a 
function of width 2a. 

We can find a lower bound to the fidelity by truncat- 
ing the sum in Ea. (|24() . By restricting L and L' to the 
domain T> = {L : ||L — R,|| 2 < A} where Rj = (x 3 \v\xj), 
we can lower bound the kernel by G(L, L') > exp{ — A^.} 
using the triangle inequality. This yields the inequality 

F( P iv,pW)>exp|-^| ^ . 

The quantity in the parenthesis is a sum over the range 
T> of a multinomial probability distribution. It is equal to 
P(||L-R|| 2 < A) > P(||L-R||i < dA) > (l- e - NdA / 2 ) 
by the Cauchy- Schwartz inequality and the typical se- 
quence theorem Eq.([Q). Thus, we get 



(25) 



Since this bound holds for all A (which is an arbitrary 
cut-off), we can maximize the RHS of Eq.J35j) — the 
optimal value turn out to be attained when A 2 = 2 ln(l + 
2Na 2 d)/Nd — to get the tightest bound: 



F(p N ,p' N ) > 1 



1 + ln(2iVcr 2 d) 
Na 2 d ' 



Hence, the measurement accuracy a can decrease as 
fast as 1/ V^V while maintaining a constant fidelity 
F(Pn,p'n) = 1 — e between the pre- and post- 
measurement states. If a decreases less rapidly than 
1/VN, e.g. N- s for < s < 1/2, the fidelity will 
go to 1 as N grows. In particular, if a is constant, 

The fidelity between the pre- and conditional post- 
measurement state, i.e. pn and p N \ e respectively, can 
be computed using similar techniques. The computa- 
tion is illustrated in the Appendix for two-dimensional 
molecules; the general case is as conceptually straight- 
forward as it is notationally cumbersome. In this case, 
even when a ^> not all measurement outcomes 

yield a state p N \ t "close" to the initial state. 
However, with a probability which approaches unity 
as 1 — exp{— ^J\og N}, fidelity will approach unity as 
1 — F(p N , Pn\i) oc exp{— vN}. In short, as N increases, 
the measurement induces a negligible disturbance to the 
state of the ensemble except with a vanishing probability. 



C. General smoothing 

We now wish to argue that the essence of our measure- 
ment accuracy-state disturbance tradeoff applies to ar- 
bitrary smoothing function qi,(i) introduced at Ea. (|l()ll . 
provided that it is actually smooth with respect to L. Let 
us be more precise. Intrinsic to the smoothing function 
is a notion of distance on the real d-dimensional vector 
space. One can define various distance measure on this 
space, e.g. our choice of smoothing function Eq. 1|11|) 
in the previous section relied on the distance ||L, L'H 2 , in- 
duced by the L 2 -norm. The exact statement of the trade- 
off will obviously depend on the choice of distance mea- 
sure. However, the essence of the result is independent 
of this choice, as all good distance measure are equiva- 
lent on small distances. Thus, a good smoothing function 
<7l(-£) should satisfy 



\9l{1)-9u{1)\<c 



|L,L'| 



(27) 



(26) 



for sufficiently small ||L — L'||i and some positive con- 
stants c and s (Eq. I|27|l is known as Lipschitz condition). 
In general, c depends on the dimension d of the molecules. 
Therefore, the dependence of the bound Eq. l|2()|) on the 
dimension d (which may seem counterintuitive) only re- 
flects our choice of the L 2 _ norm in the smoothing func- 
tion, it is not universal. Given this assumption, we can 
derive the general result. It should be mentioned that ul- 
timately, <7l (£) depend on the details of the measurement 
procedure of the corresponding macroscopic observable 
(see for example the model of Section lyijl . However, if 
this measurement is of finite accuracy, then the smooth- 
ing function must have a certain width and should satisfy 
the above assumption. 

We see from Ea. (|24|l that fidelity between the pre- 
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and post-measurement state only depends on the deco- 
herence kernel G(L,L') = / y/qL(£)qv(£)d£. Thus, the 
procedure used in the previous section carries through 
straightforwardly. We can truncate the sum Eg . 1^2411 to 
the domain T> where ||L — R||i < A. On this domain, 
the fluctuations of the kernel are bounded by Eq. Ij27(l us- 
ing the triangle inequality. Moreover, as G(L,L) = 1 by 
the normalization condition of the smoothing function, 
we obtain 



G(L,L') > l-c(A/cr) s on V. 



The bound 



F( PN ,P' N ) > {l - c I (i - e-™''*). (28) 

follows straightforwardly from the typical sequence the- 
orem Eq.Q). Given the value of c and s, one can per- 
form an optimization with respect to A to get the tight- 
est bound. However, this depends on the details of the 
smoothing function. 

Similarly, we can derive a bound for the fidelity of the 
conditional post-measurement state when a>l/ yN. In 
this case, the multinomial distribution m(L, R) behaves 
like a Kronecker delta <5l.r on the scale a, so following 

Eq. CBJ, P(Q { e N) \p N ) w q R (£). The fidelity between the 
pre- and conditional post-measurement state is therefore 



L.L' 



Again, we can truncate the sum to the domain where 
||L - R||i < A to get the bound 



F(pn,Pn\i) > 1 



«r(<) 



(1 



,-ATA 2 /2 



which goes to one with probability one as N increases, 
as illustrated in the Appendix for a specific choice of 
smoothing function. 

Finally, the scaling a ~ 1/vN is optimal. A higher 
precision would considerably disturb the state of the 
system. This is because the multinomial distribution 
m(L, R) has a width 1/y/N. Consider the expression 
of Ea.(|24p. If the kernel has a width smaller than 
the binomial distribution, the sum, and hence the fi- 
delity F(pN,p' N ), will be roughly equal to erf(a-\/iV) « 
2o-/ViV7r for a < 1/s/N. The bound is also tight for the 
conditional post-measurement fidelity F{pm , Pn\i) as n ~ 
delity is a convex function. This can also be seen in- 
tuitively by considering the behavior of two consecutive 



measurements. Upon fine grained measurement Q L 
the variance of the outcome L is 1/y/N. However, 



(AT) 



if 



we first perform a coarse grained measurement Qi 



of width a <C 1/ V^V and then perform a fine grained 
measurement on the updated state pmi, the variance of 
the second measurement outcome will be a: performing 



the coarse grained measurement has altered its statis- 
tics. This means that the coarse grained measurement 
has appreciably disturbed the state of the sample, so 
F{pN,PN\t) is far from 1. 



D. Mixed states 

The result of the previous section hold unchanged 
when the molecules of the sample are all prepared in 
the same mixed state v = Yli=i •^il'^iXV'il- The argu- 
ment proceeds in three steps. First, we can construct a 
purification of the state v 



by appending to each molecule an ancillary system of 
dimension d with orthonormal basis Clearly, the 

reduced state of the molecule — obtained by tracing out 
the ancilla — is Tr anc in a {|0)(?! , |} = v. Second, the vec- 
tors {\xj}\i}}i y j=i,... : d form a basis for the Hilbert space 
of the pair molecule+ancilla. The type projectors 
associated to the molecule only measure the type of the 
prefix Xj, so are coarse grained version of the type pro- 
jectors associated to the pair: the disturbance they cause 
to the state of the sample can only be less than the dis- 
turbance caused by the complete type projectors. Thus, 
the bound Eq.JSBJ can be applied to F(<&at,$^ v ) where 
= \cf>M® N , and 



il)$ J v(Qf ) 



1i)d£. 



Finally, by monotonicity of the fidelity 
F (£ (p) , £ (y)) > F(p,u) for any trace preserving 
quantum operation £ - the bound applies directly 
to the pre- and post-measurement state of the sample 
of molecules by tracing out the ancillas. By similar 
considerations, all of the above conclusions can be 
extended to mixed states. 



IV. EXCHANGEABILITY 

Before proceeding with the applications of the above 
results, we wish to introduce a wider class of states - 
exchangeable states — to which our results can be ap- 
plied. The concept of exchangeability was introduced in 
the classical theory of probability by de Finetti [l^ to 
substitute the incorrect use of "unknown probabilities" . 
A probability assignment is the expression of one's sub- 
jective knowledge about the possible outcomes of an ex- 
periment. Hence, it is not a property of a physical system 
itself but, rather, a property of the agent assigning the 
probability, so it can not be unknown to him! 

There are also several good reasons to believe that 
quantum states are subjective, see for example [l6l ITtI 
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ITsl IT^ | and references therein. The state of a quan- 
tum system is a mathematical construct which allows 
one to compute probabilities for various measurements 
outcomes. 1 As a consequence of the subjective nature 
of quantum states, the concept of an unknown quantum 
state is in general an oxymoron, for essentially the same 
reasons which lead to this conclusion for classical proba- 
bility assignment. 

However, unknown quantum state turn out to be quite 
useful for the description of certain physical setting. Of 
particular interest to us is the description of a sample 
of N "molecules". Under certain circumstances — e.g 
thermal equilibrium — one can arrive at the conclusion 
that all the molecules of the sample are equivalent, so 
they should all be describe by the same state v, which 
is itself unknown. This is a very common state of af- 
fair in nuclear, atomic, or molecular physics where spec- 
tral quantities — which are formally described by macro- 
scopic observables — are measured over a large collection 
of quantum systems. In fact, in almost all physical exper- 
iment where ensemble measurements are performed, the 
components of the sample are assumed to be in the "same 
unknown state" , and the purpose of the measurement is 
to (partially) determine this state. 

To arrive at an appropriate description of the sample 
without refcring to the unknown quantum state of indi- 
vidual molecules, we must clearly state the assumption of 
the agent assigning the state. His assumption is that the 
arbitrary number of molecules are all equivalent, which 
can be formalized as follows 

1. For any permutation 7r of the N molecules, ir[p] = 
p. Such a state is called symmetric. 

2. For any positive integer M, there exists a sym- 
metric state pn+m such that p^ = TrM{pN+Ai}, 
where Tru denotes the partial trace over M 
molecules. 

A state p satisfying these two conditions is called ex- 
changeable. The quantum de Finetti representation the- 
orem 0, [2(], El asserts that any exchangeable quantum 
state of a sample of N molecules can be written as 



For example, when the POVM {Ei} is measured on the 
sample, the outcome Ej is observed with probability 

P(Ej\p N ) = Tr{E jPN } 

Tr{EjV® N }Pr(v)dv 



P(EAv® N )Pr{v)dv. (30) 



We can think of P{Ej\v® N ) as the probability of Ej given 
a value of the real parameter v, but since v is unknown, 
we average this probability over the possible values of 
v distributed according to Pr(v). However, it must be 
emphasized that it is the assumption of exchangeabil- 
ity which leads to the form of Ea. (|29|) . which in turn 
legitimizes the term "unknown state" for mathematical 
convenience. 



A. Bulk tomography 

Quantum state tomography is an experimental proce- 
dures which transforms an exchangeable state of the form 



Eq. (|29*|) to a product state p 



MN 



. According to the de 



PN = I v® N Pr{v)alv 



(29) 



Finetti representation theorem, we can equivalently say 
— and this is how tomography is conventionally formu- 
lated — that the purpose of tomography is to determine 
which is the real yet unknown state v describing the N 
molecules of the sample. 

In |6(], quantum state tomography was studied in the 
context where the molecules of the sample are measured 
individually. Here, we present how quantum state tomog- 
raphy can be performed through bulk measurements. A 
similar description was recently and independently devel- 
oped in ;8j. Let An — Ylk a O) ^ e a macroscopic observ- 
able deriving from the single-molecule observable a as 
in Eq.©. Exceptionally, we denote the eigenstates and 
eigenvalues of a with a superscript a\x^) = {x^ ) 
for later convenience. The finite accuracy measurements 
of the macroscopic observable An is defined through the 

POVM {Q[ N ' a) }. 

The conditional state of the sample after the measure- 
ment of {Q^^} with outcome is 



where v are density operators of a single molecule and 
Pr(v) is a probability distribution over the quantum 
states of a single molecule. 

The interpretation of this theorem is that it is mathe- 
matically correct to look upon v as an objective element 
of reality about which we have incomplete knowledge: 
hence we assign it some probability distribution Pr(y). 



Pr{v) 



P(Q 



(N,a) 



-dv. 



PN) 



The quantity in parenthesis is proportional to the condi- 
tional post-measurement state of the sample, given that 
it was initially in state v® N . As demonstrated in Sec- 
tion IIIII this measurement has, except with a vanish- 
ing probability, very high fidelity with the original state. 
Mathematically, this means 



1 To quote Robert Griffiths, "// probabilities are not rea l, then 
pre-probabilities [quantum states] are even less reaF Il5l . 



Q f*) v ®N 



P(Q^ a) 



U ®N\®N 
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with probability approaching unity as JV grows. There- 
fore, we get 



Pn\H« 



^ N Pr{ V )P{Q[ N > a) \v® N ) 

P(Qf' a) \ P N) 



dp. 



This is as if we had updated the probability distribution 
Pr(v) of the real yet unknown state v according to Bayes' 
ruleP(j/|x) = P(x\y)P(y)/P(x). However, this is strictly 
a mathematical identity, all we did was to apply the state 
update rule Eq.© to an exchangeable state. 

We can repeat the procedure with a different macro- 
scopic observable Bn,Cn, . . . derived from the single- 
molecule observable 6, c, . . ., which do not necessarily 
commute with each other. If the sets of observables are 
sufficiently informative — i.e. if {\%^ ){ x j \}n=a,b,c,... 
contain d 2 — 1 linearly independent elements — the up- 
dated probability distribution will converge with very 
high probability to a delta function for sufficiently 
large N, Pr(v\l [a) ,l [h \ . . .) m 5(u), so the final state 
will be PN\£{a)£(b) ~ v® N . This is because the 

functions P(Q^ N '^ \v® N ) are centered around £^ — 

«4 M Vl4 rt >> (xf\v\xf), . . .) and have a width a. The 
state v is the only one satisfying all the linear constraints 
{x^\v\x^) = for all (j, = a,b,c,. . . up to accuracy 
a. Again, this is "as if the measurements simply inform 
us of the identity of the "real" v® N , without disturbing 
it in the limit of large N . 



CLASSICALITY 



Hence, the act of observing X modifies the state assigned 
to Y , However disregarding the observed value of X for 
later probability assignments is like not measuring the 
value of X at all: 



p{y J ) = Y J P{y 3 \^)P{^)- 



(31) 



Indeed, we can interpret the observation as revealing the 
"real" value of X which was there all along: the agent 
simply didn't know about it prior to her observation. 
In this sense, X — Xi IS cL real state of affair about 
which the agent learns through the act of measurement. 
Thus, the state she assign to Y prior to her observa- 
tion of X is the mixture of the state it would have given 
the different value of X, weighted by the probability of 
X, c.f. Eq. I|31|) . This reasoning extends in an obvi- 
ous way to any sequence of events X^\ X^ 2 \ . . . , X^ n \ 
We can consider that the system follows a fixed his- 
tory x^ , Xj? , . . . , x^ of which the agent has incomplete 
knowledge. 

Quantum measurements behave quite differently. A 
quantum event corresponds to a "click" on a measure- 
ment apparatus at some instant of time t. Hence, each 
event is associated a POVM element (t&) in the 
Heisenberg picture at a given time tk- (We will hence- 
forth drop the explicit time label tk-) In general, assign- 
ing definite yet unknown outcomes to these events lead 
to incorrect predictions, e.g. the sum rule 



p(e%) = J2p(e„ 



The purpose of this section is to demonstrate that 
macroscopic observables define a classical limit for closed 
quantum systems, i.e. it does not require interaction 
with any "environment" . Observations play very differ- 
ent roles in classical and quantum theory. In the classical 
setting, we can think of measurements as unveiling an 
underlying "real" state of affairs: observations reveal in- 
formation about the state of the world without affecting 
it. On the other hand, quantum measurements disturb 
or "collapse" the state of the system. 

When states, either quantum or classical, are regarded 
as subjective judgments of the world, both of the above 
descriptions need revision. Let P(xi,yj) be the joint 
probability distribution, or state, that the agent assigns 
to the classical sequence of events X = Xi and Y = yj. 
Upon the observation X = Xi, the agent updates her 
predictions for event Y according to Bayes' rule 



P(y 3 



P(xt,y 3 ) 

P(Xi) 



where P(xi) = ■ P(xi,yj). This state generally differs 
from the pre-measurement state assigned to Y 

P(y j \x i )^P(y j ) = ^P(x i ,y j ). 



does not hold in general. This is most obvious in Young's 
double slits experiment where the events E^ correspond 
to the particle going through slit j\ = 1 or 2 and % labels 
the various positions on the detector. Reasoning involv- 
ing the particle going through a definite yet unknown slit 
lead to incorrect predictions. 

There arc however sequences of quantum events which 
behave classically, as if the observations were revealing an 
underlying reality. The typical example being when all 
the POVM elements describing the events commute. The 
consistent histories approach to quantum theory |22j | lays 
down a set of conditions under which such behavior oc- 
curs. A complete list of alternatives events Cj^ 1 = {E^'} 
at time t k defines a POVM. A history is a list of POVM 
elements H = (E^ , E^J , . . . , E^) at distinct times 
t\, ti, . . . , t n . When the initial state of the system is p, 
the probability of an history H is 

P{H\p) = !>{(££>)* . . . (23«)i p (i5«)i . . . (£]«))!} 

(32) 

following Eas. (|7l9|l . A complete family of histories is the 
set of all combination of POVM elements from the sets 
C {k) at all times, T = {C (1) , h; C (2) ,t 2 ; • ■ ■ C ( ™\ *„}■ A 
family is thus a sample space on which a probability 
distribution P(H\p) is defined. The familly is said to 
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be consistent when the sum rule approximately holds for 
P(H\p). This conditions is the simplest version of all con- 
sistency conditions but will be sufficient for our purposes. 
In this sense, consistent histories define a quasiclassical 
domain of familiar experience. 

As was observed by Halliwell 01 > histories correspond- 
ing to a sequence of finite accuracy measurement of 
macroscopic observables generate a consistent family if 
the system is a sufficiently large sample of identically 
prepared molecules, i.e. when p = u® N . In other words, 
one can measure any sequence of macroscopic observable 
on a large sample and account for the observed statistics 
with a classical theory. It should be stressed that the 
single- molecule observables a, 6, c, . . . making up the his- 
tories do not need to commute. For example, the coarse 
measurement of the magnetization of a sample of spin-i 
along the z axis followed by a measurement along the 
y axis can generate a consistent family if the sample is 
sufficiently large. 

An alternative way of building our intuition in this 
direction is to consider the commutator of any two nor- 
malized macroscopic observables. Let a and b be two 
arbitrary single-molecule observables and define c to be 
their commutator c = [a, b]. These operator can be suit- 
ably normalized so that they satisfy ||a||, ||6||, ||c|| « 1. 
The normalized macroscopic observable An is defined 
as An = jj^2k=i a {k)i an d similarly for Bn and Cn', 
hence ||j4jv||, ||Bjv||, ||(7jv|| ~ 1. A straightforward cal- 
culation shows that the commutator of the normalized 
macroscopic observables obeys 

[An,Bn] = —Cn, (33) 

which implies || [Ajy, i?/v]| ~ -h- Thus, all macroscopic 
observables commute in the limit of infinite-size sample, 
and commuting observables systematically generate con- 
sistent histories: measuring the value of one observable 
does not affect the outcome statistics of other commut- 
ing observables. Once again, the infinite-sample consid- 
erations can not be applied straightforwardly to finite 
ensembles (this is the recurrent theme of this paper!) 
In particular, Eq. (|33|l does not involve any coarseness, 
which is essential to achieve consistency in finite ensem- 
bles. 

The conclusion reached by Halliwell can be extended to 
a much wider class of initial states. The first generaliza- 
tion is straightforward: by linearity of Eq. I|32l) , such fam- 
ilies are automatically consistent for initial exchangeable 
states. Indeed, for initial state pat of the form Ea. (|29|l . 
the probability of history H reads 

P{H\p N ) = J P{H\v® N )Pr{v)dv. (34) 

Clearly, if the sum rule is satisfied for the P{H\v® N ) indi- 
vidually, it is also satisfied for their convex combination. 
This is very much in the spirit of the de Finetti represen- 
tation theorem as one can interpret the outcome of the 



macroscopic measurements as revealing partial informa- 
tion about the real quantum state i/® of the sample, of 
which we have incomplete knowledge. 

Moreover, consider an arbitrary product state of the 
sample p = V\ ® v-x ® • • • ® vn- We will construct a state 
P®JV w } lose measurement outcomes, for coarse grained 
macroscopic observables, are statistically indistinguish- 
able from those obtained from the product state p. This 
will prove that Haliwell's result applies to arbitrary prod- 
uct states. Consider the symmetrized version of p: 

n[p] = jf\^2 ^( 2 ) <g> --- V < N ) 

where the sum is over all permutations of N elements. 
The reduced state of a single molecule is 

1 N 

Tr N _ x {Jl[p}} = -Y,Vk=V. (35) 
k=l 

The states H[p] and V® N are in some sense very simi- 
lar: they are both symmetric, they yield the same re- 
duced single-molecule state V, and yield the same ex- 
pectation value of the frequency operator (F^) = 
((xi|F|a;i), (xi|F|a;i), . . .). However, they are not iden- 
tical. To illustrate this, consider a sample of N two- 
dimensional molecule in the state 

p = |xi)(xi| <g) . . . ® |xi)(xi| <g> \x 2 }(x 2 \ ® . . . ® la^X^al : 

" v ' V v ' 

N/2 N/2 

half of the molecules arc in the state \x%) while the other 
half are in state |a:2). The measurement of the frequency 
operator F^ of Eq.® yields outcome (i, |) with cer- 
tainty when the state of the system is H[p] . The aver- 
age result will also be (§, 5) when the sample is in state 
V® N , but can fluctuate around this value. However, ac- 
cording to the typical sequence theorem Eq.Q, the size 
of these fluctuations will be of order 1 / vN, so can only 
be perceived by macroscopic measurements of accuracy 
a < 1/VTV. For macroscopic observables of coarseness 
a > 1/VN, the two states U[p] and V® N will yield the 
same statistics up to order a^/~N -C 1. 

Moreover, the states p and H[p] yield exactly the same 
statistics for measurement outcome of macroscopic ob- 
servables: this follows straightforwardly from the permu- 
tation invariance of the type projectors Eq.®. We have 
thus established the chain of equality 

P{Q? ] \p) = P{Q? ] \m) « P{Q[ N) \^ N l 

so the states p and 77® w yield almost identical predictions 
when a ^s> 1/ \/~N. It follows that a sequence of finite ac- 
curacy macroscopic measurements performed on a state 
of the form p = v\ (8 v<z • • • ® Vn generate a consistent 
family of histories for sufficiently large samples. 

In fact, any separable state leads to consistency of 
macroscopic histories. Indeed, when the molecules of the 
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sample are not entangled with each other, their state can 
be written as 



v x ® ui ( 



)v^Pr{vi, V2, ■ ■ ■ , v^)dv\dv2 ■ ■ ■ dv 



N- 



Now, consider the state 



T^ N Pr{v)dv 



(36) 



where we have defined 

Pr(v) = / Pr(v u . . .,i> N )d(v- — ^ Vkjdv-i . . .dv N . 
•* k 

By linearity and using the result established above, we 
have P{H\p) ~ P(H\p) provided that the histories H 
are generated by macroscopic observable of accuracy a 3> 
1 / ^/N. Since p is exchangeable, the probabilities P(H\~p) , 
and therefore the probabilities P(H\p), approximately 
satisfy the sum rule. 

To put it simply, when macroscopic measurements are 
coarse with respect to the quantum correlation length 
scale of the system, they behave classically. Indeed, 
assume that a sample of N molecules have correlation 
length £, i.e. there is ^-molecule entanglement in the 
system. Then, all of the above construction can be ap- 
plied to the N/£ collections of £ molecules. We simply 
have to treat each block of entangled £ molecules as one 
big molecule. There is no entanglement between these 
large molecules so the previous analysis applies, as long 
as the measurement accuracy is larger than y/^/N. Thus, 
we see that only entanglement on "macroscopic" scales 
can cause quantum effects to the measurement of coarse 
grained macroscopic observable. 



VI. NMR INFORMATION PROCESSING 

Room temperature nuclear magnetic resonance (NMR) 
has been for several years a benchmark for quantum in- 
formation processing [2^. The sample contains N rss 10 20 
molecules which are to good approximation non interact- 
ing due to dynamical decoupling caused by thermal ex- 
citations. Hence, the total Hamiltonian is the sum of the 
single-molecules hamiltonian H = ^2 k hrty, it takes the 
form Eq.J2J), so it is a macroscopic observable. Initially, 
the sample is in a thermal state 



,-PH 



Pn = 
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where Z = Tr{e~ 3 } and z = Tr{e~" } are the par- 
tition functions of the sample and of a single molecule 
respectively. Each molecule contains a certain number 
of nuclei which carrie a spin, and it is these spin degree 
of freedom which are used to perform the computation. 
The various spins of a molecule can have different Larmor 
frequencies ojj, which makes it possible to address them 



individually. To do so, the sample is placed by a coiled 
wire through which a sequences of externally controlled 
radio frequency (RF) current pulses can be applied. By 
properly tuning the frequency of the RF pulse, we can 
address all the spins with the same Larmor frequency, so 
all the N molecules are addressed in parallel. Therefore, 
a sequence of pulses transforms the state of the sample 
according to 



Pn 



U 



0N 



Pn 



which preserves the tensor product structure of the den- 
sity matrix p^ — v® N , it collectively changes the state 
of individual molecules v? 

It has been known for a long time [24| that the coupling 
between the nuclear spins and the coil can considerably 
disturb the state of the sample in certain regimes through 
back-action. This noise is not fundamentally irreversible, 
it is only due to our neglecting of high order terms in 
the sample's Hamiltonian. However, since the coil is also 
used to read out the state of the sample, it must unavoid- 
ably induce extra irreversible noise, of the kind discussed 
in Section IIIII This result is puzzling because the coil 
is present throughout the computation, not only during 
the measurement phase, so should in principle disturb 
the computation. 

A simple model to study the effect of this noise was pre- 
sented in Ref. . The current in the coil can be modeled 
by a continuous quantum variable Rj = J r -j\r j)(r j\dr j 
where j labels the modes of the field in the coil. Each 
field mode Rj couples to the resonant magnetization of 
the sample — i.e. to the spins of Larmor frequency ujj 
- through its conjugate momentum Pj, [Pj,Rj] — i (P 
is the "generator of translations" for R). The coupling 
Hamiltonian takes the form H c = 7 J^j PjM* where 7 is 
some coupling constant and 



N 
k=l 



N 



1)0 



(AT) 
L 



(37) 



is the total transverse magnetization of the nuclei of 
Larmor frequency ujj. Assume for simplicity that each 
molecule contain a single spin-i nucleus which couples 
to the field mode Rq . This field mode is initially in state 
\<p) = J 4>(r)\r)dr and the sample is in state pn. After a 
time t, the joint state of the field and the sample is 



Pit) 



drdr'<t>(r)<t>*(T')QVpp N QVp 
f(L')\ 



E 

L,L' 

®\r + f(L))(r' 



2 When the sequence of pulses generates a complex transformation, 
it is practically impossible to keep tract of f, as this would require 
an exponential amount of computation. Hence, given our limited 
computational capacities, the sample should really be described 
by an exchangeable state of the form Ea. 1291 . Indeed, if we 
can assign the sample a state of the form v® N after the pulse 
sequences, it means that the quantum computation was useless 
since we are able to predict its outcome! 
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where f(L) = jtN±(2L - 1). If the field mode R is 
subsequently observed to be in state |r), the state of the 
sample is updated to 



Pn 



P{r\PN) 



where Q^ = J2l< < t ) ( r ~ f{L'))Q,ffl are coarse grained 
type measurements like those of Ea. (|10fl . The width a of 
the smoothing function 4>{r— f{L')) is the initial spread of 
the field mode Rq, divided b y the coupling strengt h and 
the interaction time, a ~ y/ ((f>\R%\(f)) — {(f)\Rr,\(j)) 2 /N^ft. 
Following the results established in Section llll CI a width 
a 3> 1/VN w 1CU 10 insures us that the measurement 
does not perturb the computation. 

Of course, the measurements achieved in the labora- 
tory are much coarser than 10 -10 . Given the results pre- 
sented in this paper, we could follow and conclude 
that the presence of the coil (or the NMR measurements 
in general) induces a negligible disturbance to the state 
of the sample. However, our analysis does not apply since 
NMR measurements are not ideal (see Section lll This 
is because the coil is not in a pure state at room tem- 
perature. Consider, for example, the initial state of the 
coil in a Gaussian mixture of Gaussian-like field modes 

Pf oc J e^^ 1 \4> q ){4> q \dq where \4> q ) oc J e \r)dr 
(the 4> q are like coherent states). After a coupling time 
t, the observation of the field mode in state |r) updates 
the state of the sample to 



Pn 



where Qr^ oc 



PN\r — 



-qPNQr 



A N) 



This is the contin- 



uous version of the general state update rule Eq. (JHJ for 
non-ideal measurements. The corresponding POVM el- 

2 



ements E r — J e 



Qr 



(AO 



have width A + cr, which 



determines the accuracy of the measurement outcomes 
following Eq.0. However, the Kraus operators A qr — 

Qr_ q have width A. Following Eq.JSJ), it is this width 
which governs the disturbance caused to the state. Thus, 
it is not the measurement coarseness A + a which ulti- 
mately determines the disturbance caused to the state, 
but the details of the measurement process. 

It is therefore necessary to have a detailed model of 
the interaction between the coil and the sample to eval- 
uate its contribution to decoherence of the state of the 
molecules. We suspect that, in actual NMR settings, 
the measurement coarseness is largely due to statistical 
(thermal) fluctuations of the type of a. However, we also 
suspect the coherent spread of the coil's wave function A 
to be much larger than 10~ 10 , since coherent manipula- 
tion of the molecules appears to be possible despite the 
coupling to the coil. These questions, however, deserve a 
separate study. 



VII. CONCLUSION 

We have demonstrated a tradeoff between measure- 
ment accuracy and state disturbance for sample of identi- 
cally prepared quantum systems. A measurement coarse- 
ness smaller than 1/\N causes a disturbance to the state 
of the system which increases as the size of the ensem- 
ble grows, which is in apparent contradiction with the 
infinite-copy result. However, a measurement coarseness 
a ^> 1/ \f~N induces a negligible disturbance to the state 
of the sample. 

Using these results, we have argued that any sequence 
of macroscopic observations behave essentially classically 
provided that there is no large-scale entanglement in the 
sample. More precisely, the measurement of macroscopic 
observables generate consistent families of histories pro- 
vided that their coarseness is larger than y/£/N where 
f is the quantum correlation length scale of the system. 
We do not know however whether this scaling is optimal 
and this question deserves a separate study. Eq. ffify 
indeed suggests that consistency does not require any 
coarse-graining as N increases. The exact measurement 
of macroscopic observables indeed disturb the state of 
the sample, but might not affect the statistics of other 
macroscopic quantities. To illustrate this point, consider 
an ensemble of N two-level molecules, all prepared in the 
state |iei). Clearly, if a single molecule of the ensemble 
gets flipped to \x2), the fidelity between the original state 
and this "noisy" state is zero. In this sense, a single bit 
flip greatly disturbs the state of the sample. However, 
this single flip would have very weak repercussions on 
the statistics of the macroscopic observations. In fact, 
only a measurement of accuracy 1/N could detect this 
disturbance. The histories generated by a sequence of 
observations which do not disturb the system are guar- 
anteed to generate a consistent family. However, this is 
not a necessary condition since some type of disturbance 
— as illustrated above — do not have perceivable effects 
on the statistics of macroscopic observables. 

An interesting question arises from the study of the re- 
lation between exchangeable states and macroscopic ob- 
servations. We have seen that applying a random per- 
mutation to the molecules in a separable state yields 
a state which is not exchangeable, but possesses simi- 
lar characteristics. We do not know what type of op- 
eration can transform a generic quantum state into an 
exchangeable one. We suspect that performing a tomo- 
graphically complete set of macroscopic measurements 
on subsets of the sample followed by a random permu- 
tation of the molecules would do the trick. Physically, 
this would mean that a collective coupling to the envi- 
ronment and a diffusion process would map any state to 
an exchangeable state. This would be very interesting 
as it would extend the reach of our classicality analysis. 
Moreover, understanding under what circumstances can 
a sample of physical systems be treated as if they were 
all in the same unknown state is important since this 
is assumed in most quantum experiments performed on 
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macroscopic samples. 

Finally, we have related our study to a NMR measure- 
ment model introduced in 0. We have extended their 
analysis to the case where the coil is not in a pure mode 
state but rather in a statistical mixture of such states, 
like a thermal state. In this case, there are two parame- 
ters describing the macroscopic measurements: the width 
a of the POVM elements describes the accuracy of the 
measurements and width A of the Kraus operators gov- 
erns the disturbance caused to the state of the sample. 
Therefore, a measurement accuracy a S> 1/y/N does not 
guaranty a negligible disturbance except when the mea- 
surement is ideal. The NMR measurement process there- 
fore deserves a detailed study. 

We thanks Charles Bennett, Carl Caves, Raymond 
Laflamme, Camille Negrevergne, and Harold Ollivier for 
stimulating discussions. This work was supported in part 
by Canada's NSERC. 



VIII. APPENDIX 
A. Single molecule post-measurement state 

We will show how to compute the post-measurement 
state of a single molecule, namely Eq. For this, 

it will be useful to alter our notation a bit. In this 
subsection only, we will consider non-normalized types: 
If L(X) denotes the normalized type of X, then its 
non-normalized type is NL(X). Thus, for this section 
only, L(A") is a (/-component vector whose jth com- 
ponent Lj equal the number of occurrences of the let- 
ter Xj in X. Adding to the notation, for the type 
L = [L\, . . . , Ld) of a iV-letter string X, we denote by 
\r Xi = (Li, . . . ,Lj — 1, . . . Ld) the type of the string 
obtained by removing one occurence of Xj from X. Of 
course, this is a well defined type only when Lj > 1. 

Given this notation, we can write 



Qi N) = 



where the = when h^ Xj is not a well defined 

type. Applying the state update rule and tracing out all 
but a single molecule, we get 



(N-l) 



(38) 



Pi\t 



Tr 



N-l 



Qi N) PN Qi N) 
P(Qi N) \PN 



? £V-l)^_ lQ £V-l) 



= Y,\xi)(xiMxj)(xj\Tr v 

ij { p (Qi \pn) 

j p (Ql \Pn) 

in the last line, we used the definition of the probabil- 
ity Eq. I|15[) and the orthogonality of the type projectors 



Eq. (0J. The ratio appearing in the last line can easily 
be computed as it involves multinomial distributions, it 
is equal to 



N 



,R 



d \L1, 



N 



1 

R j 



N ' 



Inserting this in Eq. I|39|) (and keeping in mind the dif- 
ferent definitions of Lj) yields the result Eq. I|23|) . Av- 
eraging this state over measurement outcomes L gives 



Pi = ^Rj\Xj){Xj 



(40) 



The effect of coarse-grained measurements Qi on the 
state of a single molecule can be studied by straightfor- 
ward modifications of the method outlined above. The 
results are easily predictable: while the off diagonal el- 
ements \xi)(xj\ of p[ are completely suppressed when 
(7 = 0, they only get damped by a factor proportional 
to Tr{Q^[~p Q^_~p} when the measurement is coarse. 
Since ~L~ Xi and li~ Xj are very close to each other on the 
probability simplex, this decoherence factor is close to 
unity when the smoothing function qi_,(£) is sufficiently 
wide. 



B. Conditional fidelity 

We will compute, for an ensemble of N two- 
dimensional molecule initially in state (3i\xi) + /fela^), 
the fidelity between the pre- and the conditional post- 
measurement state pn and p^\i respectively. Again, we 
will use a single positive number L to denote the type 
L = (L. 1 — L). The smoothing function is Gaussian 
q L {£) = (27TCT 2 )- 1 / 2 exp{-(L - £) 2 /2a 2 }, and we assume 
that we are in the regime where a > 1/yN. Starting 
from Eq. I|17|l . we can express the fidelity 

F(PN,PN\e) = (*jv|Pjv|^|*jv) 

= Y,V<lL(t)qL'mL,\Pi\ 2 )b(L',\(3 1 \ 2 ). 

LL' 

As N becomes large, we can appeal to the central 
limit Theorem and approximate a binomial distribution 
b(L, \Pi\ 2 ) by a Gaussian distribution of mean p, = \/3i 



and variance p(l — p)/N. Moreover, the sum over the 
discrete values of L = 0, ^, . . . , 1 can be replaced by 
an integral over the range [0,1]. These substitutions are 

accurate within relative order 0(-k ) and yield 



o Jo 



27TCT 2 



(L-t)' 2 + {L' -t) 2 



N 

2-Kp(\ — p) 



-N 



(■L-Mr+(£'-Mr 

ap(i-M) dLdL'. 
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This integral can be performed and yields a complicated 
expression in terms of ref(-) function, which is not very 
enlightening. However, under the assumption that a 3> 
l/VN, one can derive the lower bound 

8 f-(^-M) 2 l 

^ JV ' p ^ ) -571^ exp l^^} 
x ( 1 - exp {-| v S}- exp {^/ I 7 Z }) 

where standard bounds on the error function were used. 
We see that as long as 



\£-p\ < A* = 0-4/2 In — (41) 

V C(T 

with c = 5\/27r/8, fidelity will be quasl-exponentially 
close to one as TV increases ( "quasi" since the correction is 
e -VN ms t eac ] f strictly e~ N ). Hence, Eq. (|41|> indicates 
for which measurements outcomes I will the conditional 
post-measurement state PN\t be quasi-exponentially close 
to p N . 

We can therefore compute the probability that the ob- 
served I falls in the bound Eq. 1)41(1 . or in other words, 
the probability that the measurement induces a negligi- 
ble error to the state of the sample. The probability of 
the measurement outcome is given by 

P(Q[ N) \ PN ) = Y J ^) h ^) ( 42 ) 

L 



following Eq. I|16(l . Using the same approximations as 
above, this gives, up to terms of relative order 0(1/ N 2 ), 

f 1 I 1 (L-i) 2 I N »(W) 2 

/ \ — /- — -e ^ci-o dL. 

Jo V 2W 2 V 27T/i(l - p) 

Finally, we can substitute the integral over [0, 1] by an 
integral over the real line; as the integrant is positive, 
this will yield an upper bound to P(Q[ N ^ \pn), which 
turns out to be a Gaussian distribution of mean p and 
variance a 2 + p(l — p)/N . Given this upper bound, we 
can ask what is the probability of getting a measurement 
outcome £ which is far from the average outcome p; 

P(\£-p\>A) = 1- T + P(Q { e N) \p N )d£ 

Jfi-A 



\^2(a 2 + p(l-p)/N) j 

__\/5A_ 

< g V5^<r 

where the last line follows from a > l/>/~N and stan- 
dard bounds on erf(-). Inserting A* of Eq. I|41|l into this 
bound gives the probability that the measurement dis- 
turbs the state of the ensemble by more than an expo- 
nentially small amount. Since a can decrease as fast as 
1 / yN, the probability of appreciably disturbing the state 
of the ensemble thus scales roughly as exp{— y / logTV}, 
and goes to zero as N grows. 
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